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The effect of pairing on small and large amplitude dynamics is discussed. Pairing correlation 
is treated in a fully microscopic transport theory using a simplified BCS version of the Time- 
Dependent Hartree-Fock Bogolyubov (TDHFB) theory. The approach is applied to study the Gi- 
ant Dipole and Giant Quadrupole Resonance in spherical and/or deformed magnesium isotopes 
showing very good agreement with QRPA theory. The present framework also provides a con- 
sistent approach to describe nuclear reactions around the Coulomb barrier. As an illustration, the 
influence of pairing correlations on multi-nucleon transfer is studied for reactions at sub-barrier 
energy. 
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1. Introduction 

Pairing correlation is known to affect significantly ground state and excited state spectroscopy 
of nuclei. On the theory side, static properties of nuclei are commonly treated using the Hartree- 
Fock-Bogoliubov (HFB) approach within the Energy Density Functional (EDF) theoiy. The Quasi - 
Particle RPA (QRPA) extend the HFB method and provides a consistent approach allowing to 
describe both low lying excited states and giant resonances. Due to the underlying small amplitude 
approximation, it cannot treat large amplitude collective motion. Large effort is currently made to 
develop the Time-Dependent HFB (TDHFB) from which the QRPA can be derived. At present, 
very few applications of TDHFB exist so f ar [jl], ^J, ||] and most of them have been made to study 
giant resonances where the QRPA already provides a competitive approach. The TDHFB theory 
is very demanding numerically and increases the computer time by a factor from 100 to 1000 
compared to the corresponding theory without pairing. This might explain why such approach has 
been up to now mainly restricted to the description of a single nucleus presenting small oscillations 
around equilibrium. In a series of workjj), || |6|], it has been shown that the BCS version of TDHFB 
can be a good compromise between the required numerical effort and the description of pairing 
beyond the independent particle picture both for giant resonances and nuclear collisions. Some 
illustrations of recent achievements are shown below. 

2. Time-dependent dynamic with pairing 

In the present approach, the many-body trial wave-function is written as a quasi-particle state 
in the BCS form 

W = II U(t)+v k (t)4(t)al k (t)) I-), (2.1) 

where u k (t) and v k (t) are the components of the special Bogoliubov transformation linking the 
quasi-particle creation/annihilation to the particle creation/annihilation of the canonical states, de- 
noted by {a k (t),a- k (t)}. This operators are associated to components (p k (r,o) in r-space and spin 
space with a\ = £ CT f drq> k (r, a)*P]j (r). Starting from the TDHFB equation and neglecting the off- 
diagonal part of the pairing field, leads to the TDHF+BCS approximation. The equations of motion 
are then conveniently written in terms of the occupation numbers n k (t) =v 2 k (t) of single-particle 
states and anomalous density components K k (t) = u* k {t)v k {t): 

iti-jj- =A* k K k - A k k£ , ifi -£■ = K t (e k - e k ) + A k {2n k - 1 ) (2.2) 

where A k is the pairing field. The advantage of TDHF+BCS compared to TDHFB is that canonical 
states evolves according to a TDHF like equation of motion: 

ih—\(p k } = (h[p}-e k (t))\(p k ) (2.3) 

where h[p] is the self-consistent mean-field while e k is a factor that is conveniently chosen as e k = 
(<Pk\h\p]\q>k). Properties of TDHF+BCS, also called Canonical basis TDHFB (CbTDHFB) have 
been extensively discussed in ref. [Q, |5|]. Such simplified approach has clear numerical advantages 
but might also lead to inconsistencies related to continuity equations ph. For this reason, it might 
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be also interesting to consider an extra simplification where occupation numbers and correlations 
are frozen in time (the Frozen Occupation Approximation (FOA)). In the following, three theories 
are compared: the TDHF approach, the TDHF+BCS and the TDHF+BCS in the FOA limit. The 
transport equations are solved on a 3D mesh using a Skyrme functional in the mean-field and a 
contact interaction in the pairing channel (for more details see ref. ^). 

3. Effect of pairing on collective vibrations in nuclei 

To study the collective response of nuclei, we consider an initial state |*P(?o)) = e -n ' F '*|^'o) 
where |*Po) is the quasi-particle ground state that is stationary solution of the BCS equation. The 
coefficient r\ should be small enough to insure the small amplitude hypothesis. The operator F 
is chosen such that specific collective modes are excited. In the following, we will consider the 
isovector Giant Dipole resonance (IV-GDR), and the isoscalar Giant Quadrupole Resonance (IS- 
GQR). Explicit form of the associated operators can be found in ref. [Q]. The nucleus response is 
obtained by solving the TDHF+BCS leading to time oscillations of F(t) = (*I / (?)|F| x I / (?)) where 



|*P(0) is the time-dependent quasi-particle state (|jj). An illustration, of the GDR oscillation is 
given in Figure [j] (top). The strength function, generally used in RPA and/or QRPA is then linked 
to the Fourier transform of F(t), denoted by F(E), through 

S(E) = ^lm(F(E)). (3.1) 

Such strength distribution is shown in Fig. [j] for the TDHF case (no pairing), TDHF+BCS 
and starting from a ground state with pairing but assuming that the pairing degrees of freedom 
are frozen in time (FOA). Note that all the calculations for giant resonances are done using the 
Skm* functional, with a mixed pairing interaction with parameters from ffijl. From the present 
comparison, one can conclude that the pairing correlations induces a modification of the GDR 
strength. In particular, a small shift of the main pic energy towards low energy is observed while 
the width of the giant resonance (life-time) is globally unchanged. Such a shift is indeed expected 
in approaches going beyond the pure independent particle picture. It is interesting to note that there 
is almost no difference between the cases where correlations are propagated in time or are frozen. 
This shows that the main effect of pairing on high-lying giant resonance is the initial fragmentation 
of the occupation numbers around the Fermi energy induced by pairing correlations. As we will 
see below this conclusion does not hold for low lying collective modes. 

As a second example, the TDHF+BCS response to a quadrupole perturbation in Mg is shown 
in figure ^[ Note that, this nucleus is found to be deformed with a deformation parameter fa = 0.33 
in its ground state. Therefore, the response is different depending on the axis direction of the 
perturbation. The response along the two main axis (K = 0) and (K = 2) are respectively shown on 
the left and right side of Fig. ^[ In both cases, the results are compared with the deformed QRPA of 
ref. [||. We see that the TDHF+BCS results, pic position and height, are in rather good agreement 
with the QRPA results although the QRPA is expected to be more general since it corresponds to the 
full TDHFB. Similar conclusion have been drawn in ref. [Q] using a different pairing interaction. 

It is worth mentioning that the TDHF+BCS results presented in Fig. ^includes the reorganiza- 
tion of correlations in time. The corresponding results where occupation numbers are frozen (FOA) 
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Figure 1: Top: Illustration of the dipole moment evolution as a function of time in 28 Mg, after an initial 
IV-GDR boost. Bottom: Strength function for 28 Mg. Both panels present three different theories, TDHF 
(dashed line), FOA (solid line) and TDHF+BCS (dots). The calculations have been performed using the 
Skyrme Skm* parameters and a mixed pairing interaction. 
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Figure 2: Isoscalar quadrupole responses for 34 Mg. Comparison between QRPA calculations [||] (solid line) 
and TDHF+BCS (dashed line). 



are shown in Fig. [|. Contrary to the IV-GDR that was not presenting low lying collective modes, 
the IS-GQR has a significant fraction of the strength at energy below 5 MeV. These low energy 
modes that are present in QRPA are absent in the FOA approximation. Therefore, while the effect 
of pairing on the high-energy collective modes could be explained by the initial fragmentation of 
the occupation numbers near the Fermi energy, the origin of low energy collective modes seems 
more complex and could only be understood through the propagation of correlation in time. 

The BCS wave function use the powerful technique of symmetry breaking, here the U(l) 
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Figure 3: Isoscalar quadrupole responses for 34 Mg. Comparison between FOA (solid line), the TDHF+BCS 
(dashed line) and the TDHF+BCS dynamics (filled circles) projected on good number of protons and neu- 
trons. 

symmetry associated to particle number conservation, to grasp pairing correlation. Accordingly, 
the response presented in Figs. [jj|3] do correspond to a weighted average of the response of nuclei 
with various proton and neutron numbers. One might worry in that case about the possible pollution 
of the strength by collective excitation of nuclei surrounding the nucleus of interest. To quantify the 
effect of particle number symmetry breaking, we have developed a projection after variation (PAV) 
method for nuclear dynamics. In this method, the response of the nucleus with proton number 
N and neutron number Z is deduced by estimating the operator F on the projected component of 

|VF(0)> i- e - 

(V(t)\P(N)P(Z)PP(N)P(Z)\V(t)) 



FNAt) (V(t)\P(N)P(Z)\V(t)) 

where P(N) and P(Z) are projectors on good particle number defined through: 

P(X) = — [ 2 *eW*-*)d(p 

2% Jo 



(3.2) 



(3.3) 



with X = N or Z. The quasi-particle states here, follow the TDHF+BCS equation of motion. 
Similarly to the original case, the strength function is obtained by performing the Fourier transform 
of Ftjpit), Such a strength is shown in Fig. ||]by filled circles. Our conclusion is that the effect of 
particle non-conservation on the collective response seems to be rather weak in general. 



4. Effect of pairing on two-nucleon transfer at sub-barrier energy 

The TDHF+BCS method can be easily implemented to perform nuclear collision using similar 
technique as for the TDHF approach. An illustration has been given in ref. [§] where the effect of 
pairing on single- and multi-nucleon transfer was analyzed. Recently, there is a renewal of interest 
regarding the transfer in the sub-barrier regime in particular to understand its possible effect as a 
competing or contributing effect to fusion [P, 10]. The understanding of pair transfer in that case is 



of special interest since it gives some information on how a composite pair of interacting particles 
might be transferred by tunneling. The present approach, by including pairing correlation as well 
as single-particle quantum dynamics provides a way to describe such process. 
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In Fig. ^ illustration of density profiles evolution during the central collision 48 Ca+ 40 Ca 
are shown for different time of the reaction. The beam energy is 48.6 MeV, that corresponds to a 
center of mass energy just below the Coulomb barrier Vg = 5 1 .9 MeV. During the collision, the two 
nuclei approaches each other up to the contact and re-separate. During the contact time, nucleons 
are exchanged as can be seen on middle and right side of Fig. |] In the TDHF framework, multi- 
nucleon exchange is treated as an independent particle process while when pairing is included, 
we do expect in an enhancement of cooperative pair transfer. This aspect could be approximately 
treated with TDHF+BCS. However, due to the break-down of the continuity equation [||], that is 
especially important when it is necessary to cut a system in two pieces, which is obviously the case 
for transfer reaction, it is more convenient to use the FOA approximation that partially cure this 
problem. The results presented below are considering this limit. Therefore, possible effects that 
are due to internal reorganization of correlations during the time evolution are neglected. 
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Figure 4: Neutron density evolution for the central collision 48 Ca (left nucleus)+ 40 Ca (right nucleus) at 
different time: before (top), during (middle) and after the contact (bottom). The time evolution of the 
density of nucleons that are initially in the 48 Ca (resp. 40 Ca) is shown in the middle (resp. right) panel. On 
the left panel, we the total neutron density is shown. 



An additional difficulty appears when theories breaking the particle number are used. Indeed, 
as in the case of giant resonances, for non vanishing pairing, the initial superfluid nucleus should 
be considered as an average over nuclei with different neutron and proton numbers. To bypass this 
problem in the case of transfer reactions, we have proposed in ref. [|]] to use a double projection 
technique. Since we are interested in the final number of nucleons in one of the side of the contact 



plane, following ref. [11 1, an operator N® = £ CT JdrW^ir, a)*I / (r, a)@(r) counting the number of 
particle in this side is defined. 0(r) is the Heaviside function equal to 1 in the considered subspace 
and otherwise. Then, one could associate to this operator, a projector 



P&(N) 



2k Jo 



2k 



Jf(N -N) 



d(p, 



(4.1) 



that gives access to the probability to find N particles in the considered space. However, contrary 
to the standard TDHF case, since the initial state does not exactly preserves particle number, one 



Effect of pairing on nuclear dynamics 



Guillaume Scamps 



should in addition define the projector onto the good particle in the total space. Accordingly, the 
probability to find N' particles in one side of the reactions while to total number of particles is 
assumed to be N is given by 



Pe(N' 



(n)\Pe(N')P(N)\n)) 

(w( t )\p(N)mt)) 



(4.2) 



The numerical details of the calculation are presented in Ref [§j. An illustration of the enhancement 
of pair transfer induced by pairing correlations is given in Fig. ||. In this figure, the one- and two- 
particle transfer probabilities obtained using TDHF and TDHF+BCS (FOA) are compared for the 
reaction 48 Ca+ 40 Ca for different energies below the Coulomb barrier. From this figure, two effects 
of pairing can be establish. First, the one-neutron transfer probability is slightly enhanced, this 
effect is due to the fragmentation of the occupation number that allows occupation above the fermi 
energy. A second effect of pairing stems directly from the non-zero initial two-body correlations. 
Those correlations increase significantly the transfer of pair of neutrons. 
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Figure 5: Illustration of one-neutron (solid lines) and two-neutron (dashed line) transfer probabilities ob- 
tained in the reaction 48 Ca+ 40 Ca as a function of the center of mass energy. The calculation are done with 
(lines) and without (lines with crosses) pairing. 



To further quantify the effect of pairing on the two-particle transfer, a systematic study of 
reactions x Ca+ 40 Ca for X = 40 to 50 has been made in ref. [||] . An illustration of the connection 
between the pair transfer enhancement and the initial pairing correlation is given in Fig. ||. In 
this figure, the ratio P2„(BCS)/P2 n {MF) for the considered reactions as a function of X and is 
systematically considered at a center of mass energy at by 6 MeV smaller than the Coulomb barrier. 
This ratio represents the enhancement of pair transfer probability due to pairing. In the top panel of 
this figure, the initial pairing gap in the x Ca collision partner is also presented. We see that the pair 
transfer enhancement is quantitatively and qualitatively strongly correlated to the initial strength of 
the pairing correlation of the considered nucleus. 
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Figure 6: Panel (a): Mean neutron pairing gap as a function of the mass for the Ca isotopic chain. Panel (b): 
Enhancement of pair transfer due to pairing as a function of the mass for the reactions x Ca+ 40 Ca. 

5. Conclusion 

In the present proceedings, several applications of the TDHF+BCS theory to collective motion 
and giant resonances illustrates how pairing can affect the small and large amplitude dynamics in 
nuclei. It is shown that pairing not only significantly affects low lying collective states as antic- 
ipated but also slightly shift the energy of giant resonances. The TDHF+BCS approach is also 
applied to study the effect of pairing on two-particle transfer. A strong enhancement that signifi- 
cantly depend on the beam energy is observed. 
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